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ABSTRACT 

The redshift of the galaxy lensing HE 0435—1223 is 0.4546 ± 0.0002, based on observations obtained 
with the Low Dispersion Survey Spectrograph 2 (LDSS2) on the Magellan Consortium's 6.5 m Clay 
telescope. Hubble Space Telescope/ ACS observations of the system also reveal a spiral-rich group of 
10 galaxies within 40" of the elliptical lensing galaxy. The redshifts for two of these galaxies were 
measured to be in the foreground (at z = 0.419) with respect to the lens, thus at least some of the 
nearby galaxies are not part of the same physical group as the lensing galaxy. Mass models of the 
system (assuming same-plane deflectors) that take the local group environment into account do better 
at explaining the observed emission-line flux ratios (which are presumably unaffected by microlensing) 
than single halo models, but the match is still not perfect. In particular, component A (a minimum 
of the light travel time) is observed to be 0.20 mag brighter than predicted and component C (also 
a minimum image) is observed to be 0.16 mag fainter than predicted. Mass models for the system 
predict an A-D time delay of either 15.8 or 17.6 days (H Q ~ 72 km s" 1 Mpc^ 1 ) depending on the 
details of the local galaxy environment. 

Subject headings: gravitational lensing: individual (HE 0435—1223) 



1. INTRODUCTION 

For a variety of reasons, follow-up observations of grav- 
itationally lensed quasars have not kept pace with the 
lens discovery rate. For example, both lens and source 
redshifts are required to convert a measured time-delay 
into an estimate of the Hubble constant (Refsdal 1964) or 
to measure the mass-to- light evolution of lensing galaxies 
(Rusin et al. 2003), but only half of the ~80 known sys- 
tems have complete redshift information. Poor knowl- 
edge of the lens galaxy environment also hinders lens- 
ing applications. Nearby groups or clusters can bias es- 
timates of the Hubble constant by contributing some 
fraction of the lensing mass convergence (Gorenstein, 
Shapiro & Falco 1988; Saha 2000), and such structures 
are likely responsible for the large (<~ 10%) shears re- 
quired to model quadruple systems (Holder & Schechter 
2003) and for explaining the observed quad-to-double ra- 
tio (Keeton & Zabludoff 2004). It is therefore important 
to understand the characteristics of each lens system as 
thoroughly as possible. In this paper, we report our mea- 
surement of the lens redshift for the quadruple quasar 
HE 0435—1223 and explore the system's local galaxy en- 
vironment from HST/ ACS imaging. 

HE 0435-1223 was discovered to be a z = 1.689 
quasar during the Hamburg/ESO (HES) survey 
for bright quasars in the Southern Hemisphere 

1 Based on observations obtained with the Magellan Con- 
sortium's Clay Telescope and the NASA/ESA Hubble Space 
Telescope. HST observations are obtained at the Space Telescope 
Science Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under NASA contract NAS 
5-26555. These observations are associated with HST program 
9744. 

2 Department of Astronomy, Ohio State University, 
Columbus, OH 43204; nmorgan@astronomy.ohio-state.edu, 
ckochanek@astronomy.ohio-state.edu 

3 Infrared Processing and Analysis Center, M/S 100-22, Califor- 
nia Institute of Technology, Jet Propulsion Laboratory, Pasadena, 
CA 91125; olga@ipac.caltech.edu 

4 Department of Physics, Massachusetts Institute of Technology, 
Cambridge, MA 02139; schech@space.mit.edu 



(Wisotzki et al. 2000). The quasar was later found 
to be gravitationally lensed by Wisotzki et al. (2002) 
during snapshot followups of HES quasars with the 
6.5 m Baade telescope at the Las Campanas Obser- 
vatory (LCO). The image morphology is a symmetric 
quad with image separations of 2'.'3 and 2'.'6 along the 
long axes of the configuration, somewhat reminiscent of 
the Northern Hemisphere "Einstein Cross" Q 2237+0305 
(Huchra et al. 1985). The lensing galaxy was promi- 
nently detected in the Baade discovery images, with 
its light profile and gri aperture colors suggestive of a 
0.3 < z < 0.5 elliptical galaxy. 

The system holds high promise for measuring a differ- 
ential time delay and estimating the Hubble constant. 
Wisotzki et al. (2002) cite evidence for both short-term 
(20% over ^2 months) and long-term (1 magnitude over 
~ 12 years) variability in the total quasar flux. The sym- 
metric image configuration also implies relatively short 
differential time delays (on the order of several weeks), 
which ought to make it straightforward to distinguish the 
time-delay signature from longer-timescale microlensing 
variability. 

Wisotzki et al. (2003) provided a refined galaxy red- 
shift estimate of z = 0.44 ± 0.02 based on integral field 
observations of the system. However, their estimate was 
obtained from low-resolution (AA = 300A) spectral re- 
binning that traced the galaxy's overall spectral energy 
distribution rather than from the detection of stellar 
absorption features. In we report a precise spec- 
troscopic measurement of the lensing galaxy redshift of 
z = 0.4546 ± 0.0002 using the Low-Dispersion Survey 
Spectrograph 2 (LDSS2) on the Clay 6.5 m telescope. 
In we make use recent HST/ ACS observations of 
the system to explore the HE 0435—1223 galaxy envi- 
ronment and find evidence for a spiral-rich group of at 
least 10 galaxies within 40", at least two of which are in 
the foreground (at z = 0.419) with respect to the lens. 
Overall the data suggest a complex lensing environment. 
We explore several lens models for the system using the 
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Fig. 1. — i2-band image of HE0435— 1223 and field obtained with the du Pont 2.5-m telescope at Las Campanas. The panel spans 7.1' 
across. Targets observed with LDSS-2 are labeled (SI, S2, lens, Gl-17). Inset: Close-up of HE 0435-1223 from Figure 2 of Wisotzki et al. 
(2002). The parallel lines show the slit orientation and width used for LDSS2 spectroscopy. 



HST/ ACS astromctry and emission-line flux ratios from 
Wisotzki et al. (2003) in 21 an d discuss the model im- 
plications and time-delay predictions in §5\ 

2. MAGELLAN /LDSS2 OBSERVATIONS 

HE 0435-1223 was observed on 2002 December f f us- 
ing the Low Dispersion Survey Spectrograph 2 (LDSS2; 
Allington- Smith et al. f994) at the Magellan Consor- 
tium's Clay telescope at LCO. The LDSS2 is a multi-slit 
spectrograph with a 7' diameter field of view and 0'.'378 
pixel -1 CCD detector. We used the medium-red grism 
blazed at 6000 A, which provided a nominal dispersion of 
5.3 A pixel -1 and a useful wavelength coverage of 4500 A 
to 8500 A. A slit width of 0'.'74 was used for all targets 
when constructing the aperture mask. The night was 
photometric and the Shack-Hartmann wavefront sensor 
helped to deliver an image quality of ^0'.'5 FWHM. 

In addition to the lensing galaxy, we obtained simulta- 
neous spectra of 19 objects within 3' of HE 0435—1223 
(see Figure 1). Our observing sequence consisted of two 
30 minute exposures with slits centered on the lensing 
galaxy and surrounding field galaxies, followed by two 10 
minute exposures through the same mask but offset by 
1'.'25 to the North-East. These last two observations were 
used to obtain quasar template spectra of components A 
and B (Figure 1 inset), which were subsequently used to 
remove quasar spillover from the lens galaxy spectrum. 
No spectrophotometric standards were observed. 



Initial data reduction consisted of bias-subtraction, 
flat-fielding, and removing cosmic-rays by interpolating 
from neighboring pixels. The two-dimensional spectra 
were then averaged together to produce a single on- 
galaxy and off-galaxy image. 

Since these are two-dimensional spectra, care must be 
taken when deriving a wavelength solution and perform- 
ing the extraction. The wavelength calibration was ob- 
tained using night sky lines bracketing each target spec- 
trum. We identified eight lines from the spectral atlas of 
Osterbrock et al. (1996), providing roughly uniform cov- 
erage between [O I]A5577 and OHA8827. The lines were 
identified on 2-4 dispersion rows (depending on the target 
slit length) with two rows always straddling the target 
spectrum, providing 16-32 points of known wavelength as 
a function of CCD row and column. We then fit a fourth 
order polynomial along the dispersion direction and a 
linear fit along the spatial direction, which provided the 
wavelength solution as a function of pixel position on the 
chip. The wavelength solution with respect to the refer- 
ence lines was extremely good: rms better than 0.8 A 
(0.15 pixels) for all spectra. The sky background was 
then removed by fitting and subtracting a line to the sky 
level for each CCD column. 

The spectral dispersion was mostly parallel to the CCD 
rows, but did bend by several pixels in the vertical from 
blue to red. To define the extraction path for the lensing 
galaxy, we fit two overlapping Gaussian profiles to each 
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TABLE 1 

Redshift Analysis for HE 0435-1223 Lens and Field Galaxies 



Object 


R.A. 


Dec. 




z 


Line Identification 


(jl 


4 38 


26.03 


-12 


18 4.6 


0.8124 


± 0.0001 


TTl /TJ a l\r\ TTT1 

[O llJ/H/3/[0 111] 


G2 


4 38 


25.25 


-12 


18 1.2 


0.3183 


± 0.0003 


r /~. TTl /it n t\/~ \ TTT1 

[O II /H/3/TO III] 


G3 


4 38 


24.85 


-12 


17 58.1 


0.3021 


± 0.0002 


\r\ TTl l\r\ TTT1 /tt 

[O llJ/[0 lllJ/Ha 


G4 


4 38 


23.54 


-12 


17 19.0 


0.3915 


± 0.0005 


\r\ TTl /tt o l\r\ TTT1 

[O llJ/H/3/[0 111] 


G5 


4 38 


22.69 


— 12 


18 57.4 


0.6240 


± 0.0007 


r /~\ TTl / X TTT1 

[O II]/[0 III] 


G6 


4 38 


22.48 


-12 


17 9.9 


0.6676 


± 0.0004 


[O II]/[0 III] 


G7 


4 38 


20.46 


— 12 


18 54.1 


0.5580 


± 0.0001 


[O IIJ/H/3 


G8 


4 38 


20.32 


-12 


16 50.7 


0.1841 


± 0.0002 


H/3/fO III1 


G9 


4 38 


18.28 


-12 


16 56.2 


0.3380 




[O II]? 


G10. . . . 


4 38 


17.73 


-12 


16 46.9 


0.3691 




[O II]? 


Gil. . . . 


4 38 


14.50 


-12 


18 29.2 


0.4579 


± 0.0001 


[O II]/H/3/[0 III] 


G12. . . . 


4 38 


24.85 


-12 


17 58.1 


0.4191 


± 0.0002 


[O II]/H/3/[0 III] 


G13.... 


4 38 


12.88 


-12 


17 2.6 


0.4189 




[O II]? 


G14. . . . 


4 38 


11.06 


-12 


16 37.7 


0.4872 


± 0.0005 


[O II] /Call H&K 


G15.... 


4 38 


10.26 


-12 


17 39.4 


0.3766 


± 0.0001 


[O H]/[0 III] 


G16.... 


4 38 


8.93 


-12 


17 29.9 


0.6880 


± 0.0001 


Call H&K 


G17.... 


4 38 


6.90 


-12 


16 37.5 


0.8434 




[O II]? 


lens .... 


4 38 


14.87 


-12 


17 14.8 


0.4546 


± 0.0002 


Call H&K 


SI 


4 38 


22.13 


-12 


17 8.3 







Star 


S2 


4 38 


19.52 


-12 


16 41.6 







Star 



Note. — Redshifts for galaxies observed with LDSS2. IDs correspond to labels in Figure 1. Redshift errors are the rms dispersion among 
multiple lines (when available), and do not include the Sz = 0.0002 uncertainty in the underlying wavelength calibration. All positions are 
J2000.0 coordinates. 



CCD column to model the flux from the A+B spectra. 
A second-order polynomial fit to A's center as a function 
of chip column gave a good fit to the dispersion drift (fit 
rms of 0.08 pixels), and was subsequently used as the 
extraction path for both the quasar and lensing galaxy 
spectra. The extraction was performed using a three 



pixel (l'.'l) width for quasar components A and B, and 
a two pixel (0'.'76) width for the lensing galaxy. For the 
remaining targets, we used a wide extraction width of 6 
pixels (2"3) parallel to the CCD rows. 

Figure 2 shows the extracted lens galaxy spectrum af- 
ter subtracting 0.2 times component A's spectrum, with 




6000 6500 7000 

Observed Wavelength (A) 



(G12, G13) 



(G1 1 , Lens) 



0.6 0.8 
Redshift 



Fig. 2. — Spectra of the lensing galaxy (bottom) and neighboring 
galaxies G12 an G13 (top and middle). The lens spectrum is shown 
after subtracting 0.2 times component A's spectrum. The four 
vertical strips denote skylines used for wavelength calibration. The 
lens, G13, and G12 spectra have been shifted by 100, 1500, and 
2250 counts, respectively, for clarity. Inset: Slice along the spatial 
direction for the lens galaxy spectrum. The two peaks correspond 
to spillover flux from the A+D and B+C quasar images, with the 
lens galaxy in the middle. 



Fig. 3. — Redshift histogram of galaxies observed with the 
LDSS2. 



4 



Morgan et al. 




Fig. 4. — PyDrizzled HST/ ACS F814W image of HE 0435—1223 and field. Prominent nearby galaxies are labeled. The concentric circles 
centered on the lensing galaxy mark 5" and 10" radii. Inset shows the inner 10" region (centered on the lens) at high contrast. 



the scale factor chosen such that the blueward side of the 
Mg II broad emission line at ~7500 A appeared smoothly- 
subtracted by eye. The Ca II H&K absorption features 
are visible, as well as the 4000 A break and G band ab- 
sorption line. Fortunately, the Ca II H&K feature fell 
150A redward of the [O I] A5577 skyline, which provided 
points for wavelength calibration on either side of the 
feature. 

To determine the galaxy redshift, we simultaneously fit 
for the the continuum plus two overlapping Gaussian pro- 
files to the Ca II H&K absorption feature and bracketing 
continuum using an appropriately fixed rest-wavelength 
separation. This does a good job of modeling the ab- 
sorption feature (reduced chi-squared of 0.7) and yields 
a best-fit redshift of z — 0.4546. The uncertainty in the 
fit is of the same order as the error in the wavelength 
calibration, ~ 1 A, translating to a redshift uncertainty 
of 0.0002. 

Figure 2 also shows the extracted spectra for galax- 
ies G12 and G13, the two closest galaxies from the lens 
that we obtained spectra for (A0=8"9 and A0=31"4 
from the lens, respectively). Both show [O II]A3727 
emission at ^5300A, and G12 also shows weak H/3 and 
[O III]AA4959,5007 emission lines at -7000A. Gaussian 
fits to the emission lines give redshifts of z = 0.4191 ± 
0.0002 and z = 0.4189 for G12 and G13, respectively. 

Table 1 lists the spectroscopic redshifts obtained for all 
18 galaxies observed with LDSS2, and Figure 3 shows the 
redshift histogram. There is no obvious group or cluster 



along the line of sight from the limited redshift sample, 
although the coincident G12,G13 redshifts at z = 0.419 
suggests some overdensity at that redshift. The red- 
shift difference between the G12,G13 pair and the lens- 
ing galaxy is 7,500 km s _1 in the G12,G13 restframe, an 
order of magnitude larger than expected for cluster asso- 
ciation and even more so for a group, so the galaxies are 
not physically associated. The closest galaxy in redshift 
space to the lens is Gil (z = 0.4579) at a projected sep- 
aration of A6>=74'.'6 (proper distance of 420 kpc), which 
is too far to significantly affect the lensing potential. 

3. HST/ ACS OBSERVATIONS 

HE 0435—1223 was observed with the Advanced Cam- 
era for Surveys/ Wide Field Camera (ACS/WFC; Ford 
et al. 1998) onboard the Hubble Space Telescope on 2003 
August 18 as part of the CASTLES imaging program 
of gravitationally lensed quasars (principal investigator, 
C. Kochanek; PID 9744). Five images each through the 
F555W and F814W filters were obtained (hereafter V 
and /) with respective total integration times of approx- 
imately 34 and 24 minutes. The data were reduced us- 
ing the IRAF/CALACS 5 package as part of the "on-the- 
fly" reprocessing at the time of download. Subsequent 
cosmic-ray rejection, geometric correction, and image 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 5.— (a): PyDrizzled H5T/F555W observations of HE 0435-1223. (6): PyDrizzled HST/F814W observations of HE 0435-1223. 
(c): Same as (b), but after subtracting the best four-component PSF model, (d): Same as (6), but after subtracting the best four-component 
PSF plus de Vaucouleurs model for the lensing galaxy. Contrast in all panels is from -0.02 e~ s~ 1 to 0.5 e~s~ 1 . 



combinations were performed using the standard PyRAF 
programs available for ACS data reduction. 

Figure 4 shows the geometrically-corrected /-band im- 
age after combining with PyDrizzle. The environment 
around the lens appears to contain a spiral-rich group 
of galaxies. There are 5 prominent galaxies within 10" 
(outer dashed circle), the closest of which (G22) is within 
5" and corresponds to the SW companion galaxy origi- 
nally noted by Wisotzki et al. (2002). Its image morphol- 
ogy is a barred spiral with arms that appear to thread 
back onto the central bulge. Galaxies G12 and G13 from 
321 are also labeled. Both show a bulge+disk image mor- 
phology, face-on in the case of G12 and edge-on in the 
case of G13, with several bright knots visible distinct 
from the central nucleus in both galaxies. The image 
morphologies are consistent with ongoing star formation 
as implied by the [O II]A3727 A emission- lines observed 
in the respective spectra. Galaxies G18, G19, G23, and 
G25 also suggest either face-on or edge-on spirals, while 
G20, G21, and G24 appear to be ellipticals. 

The V- and J-band closeup images of HE 0435-1223 
are shown in Figure 5a and 5b. The image configuration 
is consistent with previous ground-based images of the 
system, although the HST data do reveal multiple par- 
tial arcs tracing the system's Einstein radius in both the 
V- and /-band images. One can also see lensed knots 
at two different radii embedded in the arc emission and 
which likely arises from lensed structure of the quasar 
host galaxy. To obtain the relative image positions, we 



modeled the light distribution using PSFs generated with 
the TinyTim v6.1a software of Krist k, Hook (2003) and 
a circularly symmetric de Vaucouleurs profile (with ef- 
fective radius of l'/2; see below) convolved with the HST 
PSF for the lensing galaxy. The PSFs were generated 
taking into account the lens position on the WFC and 
mapped onto a 4x4 oversampled grid to assist with sub- 
pixel shifts. The relative positions and fluxes of the five- 
component model were then simultaneously solved for 
using a Powell (Press et al. 1992) minimization routine. 
The solutions were obtained in the two filters using both 
the individual (geometrically distorted) 'fit' images and 
the undistorted drizzled frames, and we found that the 
relative quasar positions agreed to within 0'.'002 using 
the different datasets. The galaxy center rms was 0'.'005 
(0'.'002) along the cardinal directions for the y(/)-band 
solutions, roughly an order of magnitude improvement 
over the earlier position obtained from ground-based ob- 
servations. 

The averaged relative offsets obtained from the 'fit' 
/-band solutions are reported in Table 2. To our initial 
surprise, the relative quasar positions differed by up to 6cr 
when compared to the Magellan offsets given by Wisotzki 
et al. (2002), even though both datasets quote similar 
astrometric precision. The differences, however, could be 
described by a 0.6% scale discrepancy between the two 
datasets, in the sense that the Wisotzki et al. (2002) 
measurements correspond to a larger angular size than 
the HST/ ACS measurements. After correcting for the 
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TABLE 2 

HST/ ACS ASTROMETRY AND PHOTOMETRY FOR HE 0435-1223 



Object A R.A. (") A Dec. (") V/V A I/Ia V-band /-band 



A =0 =0 =1 =1 18.545 ± 0.001 19.106 ± 0.001 

B -1.4772 ± 0.002 +0.5532 ± 0.002 0.607 ± 0.016 0.621 ± 0.003 19.082 ± 0.001 19.623 ± 0.001 

C -2.4687 ± 0.002 -0.6033 ± 0.002 0.579 ± 0.012 0.617 ± 0.003 19.149 ± 0.001 19.631 ± 0.001 

D -0.9377 ± 0.002 -1.6147 ± 0.002 0.557 ± 0.010 0.516 ± 0.003 19.196 ± 0.001 19.824 ± 0.001 

G -1.1687 ± 0.002 -0.5723 ± 0.002 ■ • • • ■ • 20.80 ± 0.10 19.83 ± 0.08 



Note. — Relative image positions and system flux ratios are from PSF fitting as described in the text. Quasar apparent magnitudes are from 0'/5 
aperture photometry plus a nominal 0.1 mag aperture correction. Galaxy magnitudes have been obtained using the best-fit /-band de Vaucoulcurs 
profile, with error bars corresponding to 10% uncertainties in r e . All magnitudes have been placed onto the STMAG system. 



scale difference, the Magellan and HST astrometry agree 
to better than la. 

The reconstructed /-band image without the four 
quasars (Figure 5c) allows us to investigate the detailed 
morphological properties of the lensing galaxy. We mod- 
eled the galaxy's /-band light distribution of the PSF- 
subtracted image using GALFIT (Peng ct al. 2002), tak- 
ing care to mask the ring of pixels tracing the quasar 
residuals and arc emission during the fit. The best-fit 
elliptical r 1 / 4 law gave an axis ratio of 0.83, a position 
angle of — 7? 4 East of North, and an effective radius r e 
of 1'.'20 (6.6 kpc at the lens redshift). Figure 5d shows 
the residuals after subtracting the best-fit profile. 

The effective surface brightness of the galaxy /z e , de- 
fined as the average flux inside the effective radius, is 
22.22 mag in /-band. To get a handle on the I^-band 
effective surface brightness, we applied the same /-band 
effective radius and normalized the de Vaucouleurs pro- 
file using the PSF-subtracted F-band image. This gave 
He = 23.19 mag in F-band. The total lens galaxy 
magnitude in the two filters then follows from m to t — 
H e — 5.0 log r e — 2.51og27r and are reported in Table 2. 

In Table 3, we list apparent magnitudes and relative 
positions for the 10 prominent galaxies labeled in Fig- 
ure 4 and for three nearby reference stars. Galaxy aper- 
ture magnitudes were computed using a 2'.'0 radius with 
error bars estimated from the Poisson noise inside each 
aperture. Stellar magnitudes were computed using a 0'.'5 
aperture radius and a nominal 0.1 mag aperture correc- 
tion. 

4. LENS MODELS AND TIME DELAY PREDICTIONS 

Predicting differential time delays between the quasar 
images requires an accurate mass model of the lensing 
galaxy and its environment. Wisotzki et al. (2002) found 
that a simple singular isothermal sphere (SIS) embed- 
ded in an external shear did a good job at modeling the 
quasar positions (rms of 1.3 mas), but failed to account 
for the quasar flux ratios. When compared to their i- 
band fluxes, the SIS+shear model (which only used po- 
sition constraints) underpredicted component A's flux by 
0.35 mag. Overall, it predicted components B+C to be 
almost a full magnitude (0.86 mag) brighter relative to 
components A+D than observed. 

Such discrepancies are common when modeling strong 
lenses, and it is generally accepted that the image fluxes 
are perturbed either by substructure along the line of 
sight (Dalai & Kochanck 2002) or by microlensing from 
stars in the lensing galaxy (Schechter & Wambsganss 



2002). These two processes are in principle distinguish- 
able from each other since they yield different predic- 
tions for extended sources. The flux ratios ought not 
to show a strong dependence on source size if significant 
substructure is present, but will show a source size effect 
for microlensing-induced perturbations. This is because 
larger source sizes will average over small-scale structure 
present in the microlensing caustic pattern, driving the 
system flux ratios toward the macromodel values. In 
practice, since the quasar broad-line region is several or- 
ders of magnitude larger than the continuum emitting 
region, then the emission-line fluxes ought to be much 
less sensitive to microlensing than the broad-band fluxes. 

The C IV and C III] emission-line fluxes for 
HE 0435-1223 were measured by Wisotzki ct al. (2003) 
using the Potsdam Multi- Aperture Spectrophotometer 
(PMAS) on the Calar-Alto 3.5 m telescope. They indeed 
found that the emission-line flux ratios were in better 
agreement with model predictions than were the broad- 
band values, but the match was still far from acceptable: 
component A was underpredicted by 0.26 mag and the 
combined B+C flux was still 0.57 mag brighter with re- 
spect to components A+D than observed. Attempts to 
force the SIS+shear model to reproduce the emission- 
line flux ratios could be achieved only with significant 
(~ 50a) deviations between the observed and model im- 
age positions. This difficulty could be a consequence of 
the simple isothermal model used for the lensing poten- 
tial and a more realistic model, perhaps taking the group 
environment into account, might be more successful. 

In this section, we explore several mass models for the 
system in an effort to reproduce the HST/ ACS astrom- 
etry and the Wisotzki et al. (2003) emission-line pho- 
tometry. Our basic mass model is a singular isothermal 
ellipsoid (SIE; Kassiola & Kovner (1993), Kormann et 
al. (1994)) with convergence k given by 



S c 2r \l + ecos2(9-8 e )J ' w 
where S is the projected surface mass density, E c = 
{c 2 /2-n:G){Dd/ DiDis) is the critical surface mass density, 
e is the ellipticity parameter related to the axis ratio q 
by q 2 = (1 — e)/(l + e), and the major axis is orientated 
along 9 e . The mass parameter b is given by 

>-*<(!)*£. < 2 » 

where a is the line-of-sight velocity dispersion of the dark 
matter halo. Distances /?;, D s , and Di s arc angular di- 
ameter distances to the lens, the source, and from the 
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TABLE 3 

HST/ ACS ASTROMETRY AND PHOTOMETRY FOR HE 0435-1223 FIELD 



Object 


A R.A. (") 


A Dec. (") 


V-band 


/-band 


G12. . . . 


-8.96 


3.66 


21.182 ± 0.018 


21.158 ± 0.009 


G13.... 


-30.18 


11.42 


22.486 ± 0.047 


21.980 ± 0.015 


G18.... 


-20.37 


29.56 


22.147 ± 0.043 


21.829 ± 0.012 


G19.... 


-27.69 


18.01 


21.517 ± 0.019 


21.230 ± 0.009 


G20.... 


-6.74 


13.64 


23.946 ± 0.155 


22.776 ± 0.033 


G21. . . . 


-7.62 


6.61 


22.564 ± 0.056 


21.819 ± 0.013 


G22. . . . 


-3.75 


-4.21 


22.253 ± 0.042 


21.260 ± 0.010 


G23.... 


-5.39 


-9.46 


24.409 ± 0.204 


22.850 ± 0.035 


G24.... 


1.01 


-8.39 


23.946 ± 0.155 


22.501 ± 0.024 


G25 . . . . 


-18.36 


-18.25 


22.134 ± 0.025 


21.910 ± 0.014 


SI 


9.48 


-7.16 


21.415 ± 0.004 


20.258 ± 0.002 


S2 


11.59 


-51.29 


20.183 ± 0.002 


19.250 ± 0.001 


S3 


-41.98 


67.73 


20.152 ± 0.002 


20.575 ± 0.002 



Note. — Galaxy magnitudes were obtained using a 2. aperture radius, with error bars corresponding to Poisson noise inside the aperture. 
Stellar magnitudes arc from 0''5 aperture photometry plus a nominal 0.1 mag aperture correction. All magnitudes have been placed onto the 
STMAG system. 



lens to source, respectively. For the SIS model (q — 1), b 
is also the system's Einstein radius. The two-dimensional 
effective lensing potential 4> can then be found from 
V 2 = 2k, and is simply (f> — br for the spherically sym- 
metric case (e.g., Narayan & Bartelmann 1999). We also 
consider environmental effects by adding a shear term 7 
of the form 

7 = -7yCos2(0-# 7 ), (3) 

where 7 is the shear strength and the sign convention is 
for # 7 to point toward (or away) from the mass respon- 
sible for the shear. 
We consider three models, defined as follows: 

4>ISx = + 07 (4) 

4>IEx = 4>SIE + <t>~f (5) 
4>IEISx = 4>SIE + 4>SIS,G22 + 4>-y (6) 

Model ISx is a singular isothermal sphere for the lensing 
galaxy plus external shear. It is the same model used by 
Wisotzki et al. (2002). Model IEx adds a second shear 
axis by using an elliptical isothermal halo plus external 
shear. Model IEISx explicitly takes the galaxy environ- 
ment into account. It adds a third shear axis by fixing 
a second isothermal halo at the position of the closest 
neighboring galaxy G22. 

For constraints, we use the HST/ ACS positions listed 
in Table 2 with 2 mas error circles for the quasar positions 
and a 5 mas error circle for the galaxy position. The lens 
galaxy position is allowed to vary, but its position is con- 
strained using the above error circle. The averaged C IV 
and C III] emission-line fluxes can be computed from Ta- 
ble 2 of Wisotzki et al. (2003) and give B/A, C/A, and 
D/A ratios of 0.76, 0.71, and 0.54, respectively. Wisotzki 
et al. (2003) quote formal errors for the emission-line 
fluxes on the order of 1%. While we consider models with 
1% flux errors, we also consider models with 10% flux 
errors to account for unmodclcd effects such as intrin- 
sic quasar variability. We also assume identical redshifts 
for all deflectors, even though the LDSS2 spectroscopy 
from has already showed that this is not strictly the 



case (Az = 0.036 for the lensing galaxy and G12,G13). 
Throughout this section, we adopt an (fl m ,ri\) = (0.3, 
0.7) cosmology, set the Hubble constant to 72 km s _1 
Mpc -1 , and quote all angular positions as degrees East 
of North. All models are minimized in the image-plane 
using the gravlens software of Keeton et al. (2001). 
Model results are summarized in Table 4. 

4.1. Single Halo: Models ISx and IEx 

Wisotzki et al. (2003) found that an isothermal halo 
plus external shear cannot simultaneously fit the Mag- 
ellan positions and emission-line flux ratios. We find 
a similar result using the HST/ ACS positions as con- 
straints. Using 10% flux errors, the ISx model does well 
to reproduce the image positions (9 rms = 2.7 mas), but 
predicts components B+C to be 0.59 mag brighter with 
respect to components A+D than observed. In particu- 
lar, component A's model flux is 0.28 mag fainter than 
observed. This is essentially the same result found by 
Wisotzki et al. (2003) using just the Magellan image po- 
sitions as constraints. Tightening the flux errors to 1% 
does not yield an improvement: 9 rms worsens to 5.4 mas 
and the image magnifications are mostly unchanged (see 
Table 4). 

The IEx model allows for two shear axes, external shear 
plus galaxy ellipticity. One has to be careful when solving 
for the shear and ellipticity terms since an approximate 
degeneracy exists between the two effects (Keeton et al. 
1997). To sift through the degeneracy, we solved for (b, 7, 
e) over a grid of initial conditions for (# 7 , 9 e ) and focused 
on the region of parameter space that gave the lowest % 2 
values. The best-fit model using 10% flux errors fit the 
quasar positions essentially exactly (9 rms — 0.2 mas), 
but again had trouble predicting the flux ratios, with 
the combined B+C flux 0.61 mag brighter with respect 
to A+D than observed. The problem again is chiefly 
component A, which is predicted to be 0.25 too faint. 
The galaxy parameters do agree with the observed char- 
acteristics: the model ellipticity e (defined as 1 minus 
the axis ratio) was 0.13 and pointed along — 11?5, which 
is close to the measured values of 0.17 and —7? 4 for the 
light distribution. 



8 



Morgan et al. 



TABLE 4 

Lens Model Results for HE 0435-1223 



ISx IEx IEISx 



Parameter 10% 1% 10% 1% 10% 1% 



X / d °t 


Zl.O/O 


olo/O 


n 1 1 a 
y. ( /4 


1 QA I A 

lo4/4 


O.o/o 


07 q in 

Of.O/O 


X 2 


7.1 


28.9 


0.0 


32.3 


0.0 


13.3 


riux 


8.3 


775 


8.6 


43.3 


5.5 


15.9 




5.9 


11.3 


1.1 


59.0 


0.2 


58.4 


9rms (mas) . 


O 1 

Z. ( 


K A 

0.4 




0. / 




Q V 
O. / 


Am A (mag) 


0.28 


0.28 


0.25 


0.05 


0.20 


0.03 


Am B (mag) 


-0.10 


-0.09 


-0.07 


-0.04 


0.01 


-0.01 


Amc (mag) 


-0.13 


-0.13 


-0.16 


-0.03 


-0.16 


0.00 


Am D (mag) 


0.08 


0.07 


0.12 


0.03 


0.05 


-0.02 


b (arcscc) . . 


1.20 


1.20 


1.20 


1.17 


1.17 


1.17 


b' (arcsec) . . 










0.18 


0.30 


7 


0.08 


0.08 


0.01 


0.14 


0.04 


0.14 


e 






0.13 


0.30 


0.22 


0.29 


9~[ 


-13.8 


-13.8 


-15.9 


-30.0 


-28.0 


-27.3 


". 






-11.5 


60.0 


-14.1 


60.2 


tb-A (days) 


5.48 


5.40 


7.18 


1.88 


8.72 


1.87 


t C -a (days) 


0.61 


0.44 


0.85 


-0.14 


0.70 


-0.18 


td-A (days) 


12.21 


11.89 


15.80 


4.25 


17.54 


4.07 


Note. — Model result 


3 for the 


isothermal sphere plus external shear 


(ISx), isothermal el 


lipsoid plus external shear 


(IEx), and isothermal cllipsoii 



plus external shear and isothermal sphere centered on G22 (IEISx). Models arc minimized using 10% and 1% emission-line flux errors. 9 rms gives 
the rms between predicted and observed quasar images. Mass parameters b and b' are for the lens and G22 galaxies, respectively. Magnitude 
differences are in the sense predicted minus observed. 



Tightening the IEx flux errors to 1% improved the 
flux predicitions — within 0.05 mag for all four com- 
ponents — but the galaxy parameters became unreason- 
able: the ellipticity was 0.29, almost twice as large than 
observed, and pointed nearly perpendicular (A9 = 70°) 
from the observed orientation. Moreover, the galaxy po- 
sition shifted by 0'.'038 (7-8cr) from its measured cen- 
ter. Overall the model can be rejected since it predicts 
halo properties in strong disagreement with the observed 
galaxy properties. Thus neither of the single halo models 
are consistent with both the observed galaxy properties 
and the system's emission-line flux ratios. 

4.2. Multiple Halos: Model IEISx 

The shear direction for the ISx and IEx models pointed 
roughly at —15°. Naively one would expect it to 
point toward the closest neighboring galaxy G22 at 216° 
(or equivalently 36°). The same difference was found 
Wisotzki et al. (2002) and led them to speculate that 
the neighboring galaxy might simply be a chance projec- 
tion. Another possibility is that G22 does have a signif- 
icant effect on the lensing potential, but that a second 
perturber NW of the lens dominates the shear contribu- 
tion. Figure 4 shows that this is plausible. Galaxies G12, 
G20, and G21 (as well as G18, G19 and G13 farther out) 
range from position angles of —20° to —70° and their ag- 
gregate mass or associated halo may dominate the shear 
measurement. 

Model IEISx takes both shear directions into account 
by using a halo centered on G22 for the SW perturber 
plus an external shear for the NW perturber. Using 10% 
flux errors, the model does just as well as the IEx model 



in predicting the image positions (0 rms = 0.2 mas) and 
also slightly improves the flux predictions. This time the 
B+C flux was only 0.40 mag brighter with respect to 
A+D than observed, component A was too faint by 0.20 
mag, and component C was too bright by 0.16 mag. The 
galaxy ellipticity was also slightly larger than observed, 
0.22 compared to 0.17, but did point within 7° of the 
measured orientation. The shear strength for the 10% 
model was 0.04 and pointed at —28°, or about 13° closer 
to the G20-G21-G12 group than found in the single halo 
IEx model. We again reject the 1% model since it gives 
unreasonable galaxy parameters (see Table 4) despite the 
improved image fluxes. 

5. DISCUSSION 

The improved flux predictions and closer alignment be- 
tween the shear direction and the NW group are argu- 
ments for favoring the multiple halo model, but the shear 
still does not point to any obvious perturber. It points 
between the faint G20 and G21 galaxies as seen from the 
lens. This is about +20° from the G20-G21-G12 flux- 
weighted center and about +30° from G12, the brightest 
of all neighboring galaxies within 40" of the lens (see 
Figure 6). We know that G12 is at essentially the same 
rcdshift as the lensing galaxy, so if the only significant 
perturbers were G12 to the NW and G22 to the SW, then 
we would expect the IEISx shear to point directly toward 
G12. It does not, which suggests that if the flux ratios 
are to be explained solely by the macromodel, then the 
lensing potential must be even more complicated than 
the double halo model considered here. 

Alternatively, the macromodel could be essentially cor- 
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Fig. 6. — Diagram showing relative galaxy positions (filled cir- 
cles). Dot sizes are proportional to the respective galaxy's /-band 
flux within a 2" aperture. The diagonal lines show the shear di- 
rections for the IEx (dotted) and IEISx (solid) lens models. The 
star marks the flux-weighted center of the G20-G21-G12 group of 
galaxies. 



rect and a separate effect could be perturbing the flux ra- 
tios. In general, there are four mechanisms that could al- 
ter the image fluxes: microlensing, quasar variability, dif- 
ferential dust extinction and substructure lensing. Pre- 
sumably microlensing is not an issue since we are working 
with the emission- line fluxes. Quasar variability is also 
not a factor since the A— C time delay is less than a day 
in all models, much too short compared to quasar vari- 
ability timescales (typically weeks to months). 

Wisotzki et al. (2002) already commented that ex- 
tinction is probably not very important for this system 
since the lensing galaxy is an elliptical and the Magel- 
lan colors of the four quasars agreed to <0.05 mag. In 
general, Falco et al. (1999) found median differential ex- 
tinctions among a subsample of optically-selected lenses 
of AE(B -V) = 0.04 mag. Only 10% of their sightlines 
through elliptical galaxies have 0.15 < AE(B — V) < 
0.25, so the ~0.2 mag of extinction needed to reconcile 
component C's relative faintness is unlikely but not im- 
possible. However, the PMAS resolved spectral obser- 
vations of the four quasars obtained by Wisotzki et al. 
(2003) show the quasar continuum slopes to be virtu- 
ally identical (L. Wisotzki, private communication). This 
would be a strong coincidence if any significant dust ex- 
tinction were present, and of course rules out significant 
differential extinction. 

The remaining explanation is substructure lensing. 
Millilcnsing by low-mass satellites (typically modeled as 
NFW profiles or truncated isothermal halos) in the lens 
galaxy can either brighten or dim image fluxes by sev- 
eral tenths of a magnitude. The effect has been mod- 
eled both for individual systems (Metcalf & Zhao 2002) 
and in a statistical sense (Dalai & Kochanek 2002), and 
typically requires the subclumps to comprise a few per- 
cent by mass of the primary galaxy's dark-matter halo. 



Although we do not consider substructure models here, 
such an effect remains the most plausible explanation of 
the HE 0435—1223 emission-line fluxes given the present 
data. 

We do note that substructure lensing by strictly 
isothermal halos, either embedded in the lensing galaxy 
or along the line of sight, cannot account for the ob- 
served flux ratios. Isothermal mass clumps always 
brighten positive-parity images (Keeton 2003), but com- 
ponents A and C, both positive-parity images, are re- 
spectively brighter and fainter than predicted. Of course 
the unlensed quasar flux is a free parameter, so one 
could dim the source such that both A and C became 
brighter than predicted, but then images B and D (both 
negative-parity images) would be ~0.2 mag brighter than 
predicted as well. Millilcnsing (and microlensing as 
well) tends to dim saddle-point images (Keeton 2003, 
Schechter & Wambsganss 2002), so the situation using 
isothermal halos becomes somewhat contrived. 

We can predict the differential time delays for the sys- 
tem using the measured lens redshift of z = 0.4546 ± 
0.0002 from i0 The A-D delay is the longest delay and 
therefore the most interesting to measure. The ISx, IEx, 
and IEISx models predict A-D delays of 12.21, 15.80, 
and 17.54 days, respectively. The elliptical halo (IEx) 
is formally preferred over the spherical halo (ISx) with 
X 2 /dof of 3.6 compared to 2.4, where the improvement 
comes wholly from the quasar and galaxy positions. The 
double halo (IEISx) is preferred over the single halo (IEx) 
with x 2 /dof of 2.4 compared to 1.9, where the improve- 
ment comes mostly from the images fluxes. The last 
two models yield 10% differences in the A— D time delay, 
which gives an estimate of the error that one can expect 
from uncertainties in the local galaxy environment. 

There are several ways to improve the time delay pre- 
diction. Redshift information for the remaining 8 of 10 
neighboring galaxies, especially for G22, would signifi- 
cantly improve the model confidence. In the least, it 
would point to which galaxies might share a common 
dark matter halo and thus guide the choice of paramet- 
ric models. The HST/ ACS images already show a partial 
Einstein ring for the system in /-band, and one would 
suspect a complete Einstein ring to be visible in the in- 
frared. The extra constraints offered by modeling the 
ring may be able to distinguish between the single and 
double halo models without recourse to the image flux 
ratios. For example, it would be interesting to see if the 
shear direction for a ring-constrained IEISx model points 
more in line with the center of the NW group than found 
here. 
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